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We estimate dissipative properties viz: shear and bulk viscosities of hadronic matter using rel¬ 
ativistic Boltzmann equation in relaxation time approximation within ambit of excluded volume 
hadron resonance gas (EHRG) model. We find that at zero baryon chemical potential the shear vis¬ 
cosity to entropy ratio (r//s) decreases with temperature while at finite baryon chemical potential this 
ratio shows same behavior as a function of temperature but reaches close to Kovtun-Son-Starinets 
(KSS) bound. Further along chemical freezout curve, ratio r]/s is almost constant apart from small 
initial monotonic rise. This observation may have some relevance to the experimental Ending that 
the differential elliptic flow of charged hadrons does not change considerably at lower center of mass 
energy. We further find that bulk viscosity to entropy density (C/s) decreases with temperature 
while this ratio has higher value at finite baryon chemical potential at higher temperature. Along 
freezout curve C/s decreases monotonically at lower center of mass energy and then saturates. 

PACS numbers: 12.38.Mh, 12.39.-x, 11.30.Rd, 11.30.Er 


I. INTRODUCTION 

In the context of heavy ion collision (HIC) experiments, the study of transport coefficients like shear and bulk 
viscosities is very important since these coefficients govern the evolution of the non-equilibrium system towards 
equilibrium state. In the course of nuclear collision the spatial anisotropy gets converted into momentum anisotropy 
of produced particles. The equilibration of this momentum anisotropy is governed by shear viscosity coefficient. Bulk 
viscosity relates the momentum flux with velocity gradient. Based on AdS/CFT duality authors in Ref.[lj showed 
that there is a lower bound to the value of shear viscosity to entropy ratio ( 77 /s) for any fluid found in nature. The 
lower bound (known as KSS bound) has been set to the value 1/47T. Indeed, the elliptic flow measurement at RHIC 
led to 77 /s very close to this lower boundThis discovery spurt out intense interest in transport properties of 
quantum chromodynamics (QCD). Unlike shear viscosity, bulk viscosity is very rarely studied transport coefficient. 
It was believed that the bulk viscosity does not play any significant role in the hydrodynamic evolution of the matter 
created in HIC experiments. At very high temperature QCD shows conformal symmetry whence the bulk viscosity to 
entropy ratio (C/s) is very small. But near QCD transition temperature (T c ), conformal symmetry breaking may be 
significant because the trace anomaly, (e — 3 P)/T 4 shows peak as observed in lattice simulations ))§, [|[. So one would 
expect to observe similar peak in C/s at T c as well. Indeed rise of bulk viscosity near T c has been observed in various 
effective modelsj5l4l3|. 

In order to distinguish signs of quark gluon plasma from that of hadronic matter, it is of fundamental importance 
to know the transport coefficients of hadronic matter. But computing these coefficients of hadronic matter is rather 
difficult task. There have been various efforts to compute them using some approximate schemes like Boltzmann 
transport equation in relaxation time approximation and Chapman-Enscog approximation (l3i - [27j . Some authors 
have used Green-Kubo formalism involving correlation functions and various effective models for hadronic interaction 
to compute transport properties for hadronic matter j28l - (3ll l33l]. Apart from these analytic computations shear and 
bulk viscosity coefficients have also been estimated using transport codes like ultra-relativistic quantum molecular 
dynamics (UrQMD) transport codesQ and parton hadron string dynamics (PHSD) codes[35|. 

In present work, we estimate shear and bulk viscosities of hadronic matter using relativistic Boltzmann equation 
in relaxation time approximation within ambit of (thermodynamically consistent) excluded volume hadron resonance 
gas (EHRG) model [36l - l43l ] . We assume that all hadrons have uniform hard core radius. Thus our system under study 
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is just a gas of hard sphere hadrons. As a good approximation we will use averaged partial relaxation time to compute 
viscosities. We work only in Boltzmann approximation i.e treat hadrons as classical Boltzmann particles which is 
rather good approximation at moderate temperature and take into account only elastic scattering processes. 

We will organize the paper as follows. In section II we will briefly describe the thermodynamically consistent 
excluded volume hadron resonance gas model. In section III we will give brief derivation of shear and bulk viscosities 
for multicomponent hadronic matter in relaxation time approximation and then simplify them for the case of gas of 
hard spheres assuming averaged partial relaxation time for each hadron species. In section IV we will present the 
results and discuss the implications of these results in the context of relativistic heavy ion collision. Finally we will 
summarize and conclude in section IV. 


II. EXCLUDED VOLUME HADRON RESONANCE GAS MODEL 

The only mathematically consistent statistical ensemble to describe a physical system where the particle number 
is not conserved, which is the case for relativistic heavy ion collision, is grand canonical ensemble. The knowledge of 
partition function (Z) of an ensemble is sufficient to obtain all the thermodynamical quantities of such a system. The 
pressure ( P ) of an ideal gas can be obtained from Z as 


Pideal = T lim 
V —>oo 


iTlZ-ideal (T 5 y ) 
V 


(1) 


where T is temperature, p chemical potential and V is volume of the system. In thermodynamically consistent 
excluded volume formulation where the short range repulsive interaction has been taken into account via van der 
Waals (VDW) correction to the volume one obtains the transcendental equation for the pressure as[36j 


P(T,p) = P deal (T,~p) (2) 

where p = p — vP(T,p) is an effective chemical potential with v = |47rr^ as the parameter corresponding to proper 
volume of the particle. In classical Boltzmann approximation this prescription is equivalent to additional factor of 
exp{ — vP/T ) to the pressure. Thus the pressure in excluded volume hadron resonance gas in Boltzmann approximation 
is 


P(T, p) = exp(-vP(T, n)/T)P idea i(T, p) (3) 

where Pideal in Boltzmann approximation can be written as 

Pideal{T,li) = ^ m l T2K 2^Y)cOsh(j;) (4) 

a 

where g a is degeneracy of a th hadron species and K 2 is modified Bessel’s function. We note that this expression for 
pressure is self consistent for given T and p and the pressure of interacting gas is always smaller than that of ideal 
gas. Extending HRG model by VDW type correction leads to well known suppression of particle number densities 
consistent with heavy ion collision experiments [!i|. Other thermodynamical quantities can be readily obtained from 
pressure, viz: Baryon number density Ub = dP(T, p)/dp, entropy density s(T,p) = dP(T,p)/dT, energy density 
e(T, p) = Ts(T , p) — P(T, p) + pub{T, p), sound speed C%(T , p) = dP(T, p)/de(T, p). 


III. DISSIPATIVE PROPERTIES IN RELAXATION TIME APPROXIMATION 

In this section we will follow method given in Ref. j44j to calculate transport properties in relaxation time approxi¬ 
mation. Fundamental equation of the kinetic theory is Boltzmann transport equation given by 
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v p = P/Ep is single particle velocity and I{f p } is collision integral and gives the rate of change of distribution function 
f p due to collisions. f p is in general non-equilibrium distribution function. 

Relaxation time approximation assumes that the collisions always bring the system to the local equilibrium expo¬ 
nentially with relaxation time which is of order of collision time. Thus in this approximation collision integral can be 
written as 


Up - f° P ) 


t(E p ) 


( 6 ) 


where t(E p ) is called relaxation time or collision time which in general depends on energy of the particle and f p is 
equilibrium distribution function given by 


f° = 

J P 


exp 


Ep—p.u—fi 


± 1 


( 7 ) 


where u is fluid velocity and ± corresponds to fermions and bosons respectively. In hydrodynamical description of 
QCD matter shear and bulk viscosities enters in dissipative part (T^T) of stress energy tensor. 


rpfLV _ rpfiV _|_ rj-iflV 


( 8 ) 


where Tq V is ideal part of stress tensor. In the local lorentz frame dissipative part of stress energy tensor can be 
written as 


± dissi 




( 9 ) 


In terms of distribution function T^ issi can be written as 


T lJ = 

± dissi 


(Syj-Wp 


( 10 ) 


where 5f p is the deviation of the distribution function from equilibrium which governs the dissipative properties of 
the system. From Eq. 10 and Eq. © we get 


dfO df 0 

5f p = -r{E p )[-±+v' p Jz. 


dt 


' dx l 


( 11 ) 


Assuming steady flow of the form u l = (u x (y), 0,0) and space-time independent temperature, Eq. © simplifies to 
T xy = — pduxjdy. From Eq. (fTUl) and Eq. (TTTT) we get (using Eq.© with p = 0) 




Thus coefficient of shear viscosity for a single component of hadronic matter is finally given by 

1 /' d3 P \P A f o 

V 15 T J (27r) 3T( p ^E? fp 


( 12 ) 


(13) 


Bulk viscosity is related to the dissipation in the system when it is uniformly compressed. Taking trace of Eq. 
we get 


C TdissiK = -3C^7 


(14) 


Also from Eq. ms and Eq. ED we get 


" 2 'df° 


;5/,° 


(T dlssi )l - J (2 ^ 3 t{E p ) P ^ + v' p 


(15) 
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Using energy momentum conservation law = 0, together with Eq. (fT4|) and Eq. (fl~5|) one can arrive at [44 


c = I 


d 3 p 




KCi - 


T J (2nf y p,j p[ p nB 3E p _ 

where C% B = ^| nfl , is speed of sound at constant baryon density. 


(16) 


Thus for multicomponent hadron gas at finite chemical potential shear and bulk viscosities can be written as 

d 3 p p 4 


V = 


-V 


15T 


(2tt)3 El 


(.Tafa + Taf°) 


(17) 


c = ^E 

a 

+ f afa 


d 3 p 

w? 


T f 

1 aj a 


+(b -- 


3 E a 


fc 2 i _zE 2 

“ nB ' dn n ) e 3 E a _ 


(18) 


where E'l = p 2 + ml and ( dP/dns)e = n_e/ [dnsfdp) + Cl B T 2 d{p,/T) / dT. In above expressions bar stands for 
contribution of antiparticles. 

Energy dependent relaxation time is defined by expression 


L/ P n / d3 Pb d3 Pc d3 Pd VtT( tv A\ fO 

( E a) = B I ^ c ’ d )fb 


bed 


(27r) 3 (27 t) 3 (27t) 3 


M 


where the transition rate W(a,b c, d) is defined by 

W(a, b -+ c, d) = Wifa+W-Pc-Pd) 

1 ’ ’ 1 2E a 2E b 2E c 2E d 

with | M | being transition amplitude. In the center of mass frame Eq. TO can be simplified as 

-ln> V" f d3 Pb - 4m2 f o _ V" f d3 Pb 

T (Ea) = 1 _^J J^CTab 2 E„ 2 Eh f»= 2 ^ I 


(2tt) : 


7&a b Vabf b 


(19) 


( 20 ) 


( 21 ) 


where v ab is relative velocity and \fS is center of mass energy. a ab is the total scattering cross section for the process, 
a{po) + b{p b ) ->• a(p c ) + b{p d ). 

For the simplicity we can use averaged relaxation time (f) which is rather a good approximation as energy dependent 
relaxation time[45}. One can obtain f as follows. Averaging over / 3 Eq. m becomes 


Sifrr-HEaW = ^jJ^^ ab v ab f° a fS 


r d3 P a f o 

J (2k) 3 J a 


(2tt) 3 -la 

Thus averaged partial relaxation time is given by 

T ” 1 = E n bWabVab) 


r d 3 p a f 0 
J (2k) 3 da 


( 22 ) 


(23) 


where n b = f j$pf b is the number density of b th hadronic species. 

In this work we will use equilibrium Maxwell-Boltzmann distribution (in the local rest frame) given by 


fa = exp[ - 


Eg - n 
T 


(24) 


The thermal average of total cross section times relative velocity i.e (av) for the scattering of hard sphere particles 
(having constant cross section, a) of the same species at zero baryon density can be calculated as follows [46L [47j | . 
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With Maxwell-Boltzmann distribution f°(E) = exp(—E/T), the thermal average ( av) for the process a(p a ) + a(p b ) —► 
a(p c ) + a(pd) can be written as 


i&abVab) — 


■/ d 3 p a d 3 p b v ab e E °-/ T e Eb/T 


f d 3 p a d 3 p b e~ Ea / T e~ Eb / T 
Momentum space volume elements d 3 p a d 3 p b can be written in terms of scattering angle 9 as 


(25) 


d 3 p a d 3 p b = (4Tr) 2 p a p b dE a dE b ^cosd (26) 

Changing integration variables from E a ,E b ,9 to E_,E + ,S we gets 

d 3 p a d 3 pb = 2tt 2 E a E b dE_dE+dS (27) 

where S = ( p a +p b ) 2 is usual Mandelstam variable and E± = E a ± E b . With this change in variables the integration 
region transform as 


E _ < (28) 
with E + ^ y/S and S ^ 4m 2 . Thus numerator in Eq. 03 becomes 

J d 3 pad 3 p b v ab e- E °-l T e~ Eb / T = 2t t 2 T J dSVS(S - Am 2 )Ki(\fS/T) (29) 

Similarly denominator of Eq. (1251) can be evaluated as 

J d 3 p a d 3 p b e~ Ea ^ T e~ Eb ^ T = [47rm 2 TJ'L 2 (m/T)] 2 (30) 

where K n is modified Bessel function of order n. Thus the thermal average becomes 

/•OO 

{°abVab) = R iT ° ( m / dSVS(S - 4,m 2 )K\(\fS/T) (31) 

8m 4 TK|(m/T) J im 2 

For scattering between different species of the particles ( a(p a ) + b{p b ) —>■ a(p c ) + b(pd)) one can generalize above 
equation to get 


{^ab^ab) — 


8Tmlm 2 b K 2 (^)K 2 (^) 


' m a -\-mb 


[5-(m a -m fe ) 2 ] 
dS - -js - [S - [m a 


m b ) 2 ]Ki(y/S/T) 


(32) 


Computing the thermal averaged cross section as above, one can relate it to the relaxation time in Eq. E51) . The 
viscosities can then be calculated using Eq. ea and Eq. m once the thermodynamic quantities are estimated using 
EHRG model. 


IV. RESULTS AND DISCUSSION 

To estimate the thermodynamical quantities using EHRG model we take all hadrons and their resonances with mass 
cut-off 2.25 GeV[48] in the partition function. The only unknown parameter needed to compute viscosity coefficients 
is hadronic hard core radius r b . On one hand in Ref. [491] it has been argued that one can take uniform hard core radius 
equal to 0.3fm for all hadrons. This argument is based on the hard core repulsive interaction between nucleons. For 
mesons there is no detailed results for short range repulsive interaction. But one can set same hard core radius to all 
mesons based on similarity of charge radii compared to baryons and pion-nucleon phase shift energy dependence. On 
the other hand in Ref. [50] authors demonstrated that strangeness Horn behavior can be described in hadron resonance 
gas model without spoiling hadron yield fit if one set vanishing hard core radius (r^) for pions, rh = 0.35fm for kaons, 
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0.3fm for all other mesons and 0.5fm for all baryons. Based on all these arguments we take two values of 77 0.3fm 
and 0.5fm for all hadrons to compute viscosity coefficients. The cross section of the hadrons in terms of hadron hard 
core radius is given by a = 47rr^. 

Results for shear viscosity coefficient is shown in Fig. ©■ The general behavior of 77 /s is similar to that ob¬ 
served in Ref. [Tsl] where the authors considered EHRG within relativistic molecular kinetic theory unlike relaxation 
time approximation scheme used in this work. We also compare our results of 77 /s with other model calculations 
like Chapman-Enscog theory [5l]], scaling hadron masses and couplings (SHMC) model [52t| and chiral perturbation 
theory [531] at zero baryon chemical potential as shown in Fig. ©■ We note that the general behavior of 77 /s is in 
conformity with these models. We also note that at low temperature (~ 0.120GeV) where the pions are dominating 
degrees of freedom, our results matches with Ref. [5^] where the authors estimates 77 /s for the gas of pions using 
chiral perturbation theory while at high temperature (above 0.120 GeV) our results matches with Ref. [52] where 
the authors estimated this ratio in SHMC model for hadronic matter. Further we observe that at finite chemical 
potential, although the general behavior of the ratio is similar as a function of temperature, ratio is smaller than that 
at 77 = OGeV and approaches closer to KSS bound. Thus finite baryon chemical potential significantly affect 77 /s. 
Although the shear viscosity itself increases with 77 as shown in Fig.©), decrease in ratio 77 /s at finite 77 is solely due 
to rapid increase in entropy density. This behavior of 77 /s at finite baryon density is consistent with Ref.(5l} where 
the authors have estimated 77 /s using Chapman-Enskog theory within hadron resonance gas model. We compare our 
results with the results of Ref. [51} at zero and finite 77 as shown in Fig.(©). We note that the general behavior of 77 /s 
is similar except the fact that the value approaches closer to KSS bound in Chapman-Enskog theory. 

It is important to note that at finite chemical potential, 77 /s cannot be inferred as a measure of fluidity [54^. Also 
this ratio can be shown to violate KSS bound in kinetic theory. Based on crude kinetic theory argument one can 
show that 77 ss A y/„(?i(p)A) a , where n is number density, (p) is thermal momentum and A is mean free path. Kinetic 
theory is valid only for those gases for which mean free path is much smaller than the typical size of the system (L) 
i.e, A < L and for those gases for which A must be larger than inter particle spacing. Then uncertainty relation 
A (p) > 1 implies that there is a lower bound to shear viscosity, 77 > 2T 3 [5|. Also in the non-relativistic limit one 
can show that the shear viscosity of the gas of hard spheres is independent of number of particle species [l4]. On the 
other hand, entropy density of the gas consisting multiple hadronic species (which also goes as T 3 ) can be made very 
large so that ratio 77 /s can be made arbitrarily small. In fact at sufficiently high chemical potential mixing entropy 
of multicomponent hadron gas overwhelms and hence ratio 77 /s can go below KSS bound. This fact has been used in 
Ref. [5^| to give counterexample to KSS bound. We might also like to mention here that KSS bound can be violated 
in certain field theories [57H59f . 

At finite 77 , it is the ratio 77 T/(e + P) which is correct measure of fluidity[54j|. Quantity e + P is called enthalpy 
and as per thermodynamical relation, e + P = Ts + 77773 , we note that at 77 = 0 we get back 77 /s as a fluidity 
measure. From Fig. ([3]) we note that effect of finite chemical potential is more pronounced in ratio r]T/(e + P). This 
can again be attributed to rapid rise in enthalpy at finite 77 . The general behavior of the ratio r\T /(e + P) is again in 
conformity with Ref.JlJ as shown if Fig. (Oo) except the fact that for given chemical potential the ratio is smaller in 
Chapman-Enskog theory. 

Fig.© shows results for the bulk viscosity. We note that ratio (/s decreases with temperature at zero chemical 
potential. As shown in Fig.©, the general behavior of C/s is similar to that observed in Ref.[52] where the authors 
estimated the bulk viscosity using SHMC model. At finite chemical potential although the ratio (/s decreases at low 
temperature, it increases in the window T = 0.120 — 0.160 GeV. This is because bulk viscosity itself increases very 
rapidly in this window as shown in Fig.©). This rise may be attributed to the explicit scale symmetry violation by 
finite chemical potential and hence the massive nucleon excitations which contribute more at finite baryon chemical 
potential[29|. We might mention here that alth oug h the inelastic scattering processes needs to be taken into account 
for the precise estimation of the bulk viscosity |60f. authors in Ref. showed that inelastic processes are irrelevant 
in the bulk viscosity computation at low and moderate temperatures. In Fig.© we compare (/s estimated in our 
model with SHMC model[52||. We note that our (/s curve vanishes faster at high temperature 

One can make connection with heavy ion collision experiments by finding the beam energy (y/S) dependence of 
the temperature and chemical potential. This is extracted from a statistical thermal model description of the particle 
yield at various yfS [62(. The freeze out curve T(fi) is parametrized byT(77) = a — bp, 2 — C77 4 , where, a = 0.166 ± 0.002 
GeV, b = 0.139±0.016 GeV -1 and c = 0.053±0.021 GeV -3 . The energy dependence of the baryon chemical potential 
is given as 77 = d/( 1 + eyfS), with, d = 1.308 ± 0.028 GeV, and e = 0.273 ± 0.008 GeV -1 . From Fig.©) we observe 
that ratio 77/s is well above KSS bound at low center of mass energy and increases monotonically to become constant 
at higher y/S along freezout curve. This is legitimate since low y/S corresponds to low temperature and high chemical 
potential along freezout curve at which shear viscosity is smaller. Fig.©>) shows ratio r]T/(e + P) along chemical 
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freezout. We observe that this ratio again remains constant apart from initial rise. Since ratio r]T/(e + P) is a true 
measure of fluidity at finite baryon chemical potential, we conclude that in chemical freezout transition the fluid 
behavior of hadron gas does not change [5l|. Further, along freezout curve ratio (,/s decreases monotonically first and 
then becomes independent at higher center of mass energies as shown in Fig. dZ}c). 


V. SUMMARY AND CONCLUSION 


In this work we have estimated dissipative properties of hot and dense hadronic matter using Boltzmann equation 
in relaxation time approximation within ambit of excluded volume hadron resonance gas model. We assumed uniform 
hard core radius to all hadrons and included all hadrons and their resonances with mass cut-off 2.25 GeV. We 
considered only elastic scattering processes to compute transport coefficients. We find that shear viscosity to entropy 
ratio decreases with temperature. Further at finite chemical potential 77 /s shows same behavior as a function of 
temperature but ratio is smaller as compared to /r = 0. This decrease is solely due to rapid increase in entropy density 
at finite /i. As pointed out in Ref. [541], at finite baryon density it is r]T/(e + P) and not 77 /s as a correct measure of 
fluidity. We find that effect of finite /i is more pronounced for rjT/(e + P) and this is again attributed to rapid rise 
in enthalpy. Further the bulk viscosity to entropy ratio decreases with temperature but at higher baryon chemical 
potential this ratio is higher as compared to zero chemical potential. This is due to heavier baryonic excitations 
which makes additional contribution to the bulk viscosity at finite baryon chemical potential. Further, along chemical 
freezout curve both the ratios 77 /s and t]T/(e + P ) remains constant apart from initial rise. This suggest that fluid 
behavior of hadron gas does not change along chemical freezout transition. Further, initially along along freezout line 
ratio C/s decreases monotonically and then becomes independent of central of mass energy . 

In this work we have taken zero temperature hadron masses independent of temperature (T) and chemical potential 
(/i) in the partition function. But hadron masses depends on T and n since constituent quark masses (which makes 
large contribution to hadron mass) arises due to interaction with chiral condensates. These chiral condensates are 
depends on T and fi and actually melts as the temperature and chemical potential increases. So it would be more 
realistic to include T and fi dependent hadron masses in HRG model and to see the thermodynamics of such system 
and its effect on transport properties of hadronic matter. Work in this direction is in progress and will be reported 
elsewhere (63j|. 
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FIG. 1: (Color online) Top panel shows shear viscosity to entropy density ratio 77 /s for = 0.3 fm as a function of temperature 
for different chemical potentials. Middle panel shows rj/s for rh = 0.5 fm. Bottom panel shows the shear viscosity as a function 
of temperature at different chemical potentials for rh = 0.5 fm. 
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FIG. 2: (Color online) Comparison of rj/s estimated in our model with the estimations of other models f5H - f53il at zero chemical 
potential. 
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FIG. 3: (Color online) Left panel shows ratio r]T/(t + P ) for rn = 0.3 fm as a function of temperature for different chemical 
potentials. Right panel shows r/T/(e + P) for rn = 0.5 fm. 

















T (GeV) T (GeV) 

(a) (b) 


FIG. 4: (Color online) Left panel shows comparison of r/Js (rh = 0.3fm) estimated in our model with the estimations of Ref. 
Right panel shows comparison of r/T /(e + P) with Refi51]. 
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FIG. 5: (Color online) Top panel shows bulk viscosity to entropy ratio (/s for rh = 0.3 fm as a function of temperature for 
different chemical potentials. Middle panel shows (,/s for rh = 0.5 fm. Bottom panel shows the bulk viscosity as a function of 
temperature at different chemical potentials for rh = 0.5fm. 
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FIG. 6: (Color online) Comparison of C/s estimated in our model with estimations using SHMC model[5^] at zero chemical 
potential. 
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FIG. 7: (Color online) Top panel shows 77 /s along freezout line for = 0.5fm. Middle panel shows r]T/(e + P) along freezout 
line. Bottom panel shows (,/s along freezout line for = 0.5fm. 






